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Abstract. We use the evolutionary turbulent model of Jupiter's subnebula described by Alibert et al. (2005a) 
to constrain the composition of ices incorporated in its regular icy satellites. We consider CO2, CO, CH4, N2, 
NH3, H2S, Ar, Kr, and Xe as the major volatile species existing in the gas-phase of the solar nebula. All these 
volatile species, except CO2 which crystallized as a pure condensate, are assumed to be trapped by H2O to 
form hydrates or clathrate hydrates in the solar nebula. Once condensed, these ices were incorporated into the 
growing planetesimals produced in the feeding zone of proto- Jupiter. Some of these solids then flowed from the 
solar nebula to the subnebula, and may have been accreted by the forming Jovian regular satellites. We show 
that ices embedded in solids entering at early epochs into the Jovian subdisk were all vaporized. This leads us 
to consider two different scenarios of regular icy satellites formation in order to estimate the composition of the 
ices they contain. In the first scenario, icy satellites were accreted from planetesimals that have been produced in 
Jupiter's feeding zone without further vaporization, whereas, in the second scenario, icy satellites were accreted 
from planetesimals produced in the Jovian subnebula. In this latter case, we study the evolution of carbon and 
nitrogen gas-phase chemistries in the Jovian subnebula and we show that the conversions of N2 to NII3, of CO to 
CO2, and of CO to CH4 were all inhibited in the major part of the subdisk. Finally, we assess the mass abundances 
of the major volatile species with respect to H2O in the interiors of the Jovian regular icy satellites. Our results 
are then compatible with the detection of CO2 on the surfaces of Callisto and Ganymede and with the presence 
of NII3 envisaged in subsurface oceans within Ganymede and Callisto. 
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1. Introduction 

' In a recent paper, Alibert et al. (2005a) (hereafter referred 
to as Paper I) developed a two-dimensional time depen- 
dent a-turbulent model of the Jovian subnebula whose 
evolution is ruled by the last sequence of Jupiter forma- 
tion. These authors carried out migration calculations in 
the Jovian subnebula in order to follow the evolution of 
the ices/rocks ratios in the protosatellites as a function of 
their migration pathways. By tempting to reproduce the 
distance distribution of the Galilean satellites, as well as 
their ices/rocks ratios, they obtained some constraints on 
the viscosity parameter of the Jovian subnebula and on 
its thermodynamical conditions. 

They showed that the Jovian subnebula evolves in two 
distinct phases during its lifetime. In the first phase, the 
subnebula is fed through its outer edge by gas and gas- 
coupled solids originating from the protoplanetary disk 
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as long as it has not been dissipated. During the major 
part of this period, temperature and pressure conditions 
in the Jovian subnebula are high enough to vaporize any 
icy planetesimal coming through. When the solar nebula 
has disappeared, the subnebula enters the second phase 
of its evolution. The mass flux at the outer edge stops, 
and the Jovian subnebula gradually empties by accreting 
its material onto the forming Jupiter. At the same time, 
due to angular momentum conservation, the subnebula 
expands outward. Such an evolution implies a rapid de- 
crease of temperature, pressure and surface density condi- 
tions over several orders of magnitude in the whole Jovian 
subnebula. 

In the present work we focus on the possibility of esti- 
mating the composition of ices incorporated in the regu- 
lar icy satellites of Jupiter in the framework of the model 
described in Paper I. A similar study was previously con- 
ducted by Mousis & Gautier (2004) (hereafter referred to 
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as MG04) but here we present several significant improve- 
ments. 

First, the initial accretion rate of our turbulent model 
of the Jovian subnebula is fully consistent with that cal- 
culated in the last phase of Jupiter formation (see Paper 
I). As a result, the temporal evolution of our model of the 
Jovian subnebula, as well as the thermodynamical condi- 
tions inside the subnebula, are quite different from that of 
MG04. Hence, the question of the resulting composition of 
ices incorporated in the Galilean satellites remains open. 

Second, in our model, the solids flowing in the sub- 
nebula from the nebula were formed in Jupiter's feeding 
zone. For the sake of consistency, it is important to cal- 
culate their composition using the same thermodynamical 
and gas-phase conditions as those considered by Alibert 
et al. (2005b - hereafter referred to as A05b). Indeed, 
using the clathrate hydrate trapping theory (Lunine & 
Stevenson 1985), A05b have interpreted the volatile en- 
richments in Jupiter's atmosphere, in a way compatible 
with internal structure models derived by Saumon and 
Guillot (2004). As a result, they determined the range of 
valid C02:CO:CH4 and N2:NH3 gas-phase ratios in the so- 
lar nebula to explain the measured enrichments, and the 
minimum H2O/H2 gas-phase ratio required to trap the 
different volatile species as clathrate hydrates or hydrates 
in icy solids produced in Jupiter's feeding zone. Our cal- 
culations then allow us to determine the composition of 
ices in Jupiter's regular satellites, in a way consistent with 
the enrichments in volatile species observed in the giant 
planet's atmosphere by the Galileo probe. 

Finally, we consider further volatile species that are 
likely to exist in the interiors of the Jovian regular icy 
satellites. In addition to CO, CH4, N2, and NII3 that 
have already been taken into account in the calculations 
of MG04, we also consider CO2, Ar, Kr, Xe and II2S. 
CO2 has been detected on the surface of Ganymede and 
Callisto (McCord et al. 1998; Hibbitts et al. 2000, 2002, 
2003) and is likely to be a major carbon compound in the 
initial gas-phase of the solar nebula since large quantities 
are observed in the ISM (Gibb et al. 2004). Moreover, Ar, 
Kr, Xe and H2S abundances have been measured in the 
atmosphere of Jupiter (Owen et al. 1999). Since, accord- 
ing to A05b, these volatile species have been trapped in 
icy planetesimals in Jupiter's feeding zone during its for- 
mation, they may also have been incorporated into the 
material (gas and solids) delivered by the solar nebula to 
the Jovian subnebula and taking part in the formation of 
the regular satellites. 

The outline of the paper is as follows. In Sect. 2, we ex- 
amine the conditions of volatiles trapping in solids formed 
in Jupiter's feeding zone. In Sect. 3, we recall some details 
of the thermodynamical characteristics of our turbulent 
model of the Jovian subnebula. This allows us to inves- 
tigate the conditions of survival of these solids formed 
inside the solar nebula and accreted by the subnebula. 
In this Section, we also study the evolution of the gas- 
phase chemistries of carbon and nitrogen volatile species 
in the subdisk. In Sect. 4, we estimate the mass ratios 



with respect to water of the considered volatile species in 
the interiors of regular icy satellites. Sect. 5 is devoted to 
discussion and summary. 

2. Trapping volatiles in planetesimals formed in 
Jupiter's feeding zone 

The volatiles ultimately incorporated in the regular icy 
satellites have been first trapped under the form of hy- 
drates, clathrate hydrates or pure condensates in Jupiter's 
feeding zone. The clathration and hydratation processes 
result from the presence of crystalline water ice at the 
time of volatiles trapping in the solar nebula. This latter 
statement is justified by current scenarios of the formation 
of the solar nebula who consider that most of ices falling 
from the presolar cloud onto the disk vaporized when en- 
tering in the early nebula. Following Chick and Cassen 
(1997), II2O ice vaporized within 30 AU in the solar neb- 
ula. With time, the decrease of temperature and pressure 
conditions allowed the water to condense and form micro- 
scopic crystalline ices (Kouchi et al. 1994, Mousis et al. 
2000). Once formed, the different ices agglomerated and 
were incorporated into the growing planetesimals. These 
planetesimals may ultimately have been part of the ma- 
terial (gas and solid) fiowing in the subnebula from the 
solar nebula. Moreover, larger planetesimals, with metric 
to kilometric dimensions, may have been captured by the 
Jovian subnebula when they came through. 

In the model we consider, Jupiter forms from an em- 
bryo initially located at ~ 9-10 AU (Ahbert et al. 2005c). 
Since the subnebula appears only during the late stage of 
Jupiter formation when the planet has nearly reached its 
present day location (see Paper I for details), we used the 
solar nebula thermodynamical conditions at 5 AU. On the 
other hand, the use of the solar nebula thermodynamical 
conditions at ~ 10 AU would not change our conclusions, 
since the composition of icy planetesimals does not vary 
significantly along the migration path of Jupiter if a sim- 
ilar gas-phase composition is assumed (A05b). 

The trapping process of volatiles, illustrated in Fig. 
n is calculated using the stability curves of clathrate hy- 
drates derived from the thermodynamical data of Lunine 
& Stevenson (1985) and the cooling curve at 5 AU taken 
from the solar nebula model used to calculated Jupiter's 
formation (see A05b). For each considered ice, the domain 
of stability is the region located below its corresponding 
stability curve. Note that the use of cooling curves derived 
from others evolutionary a-turbulcnt models of the solar 
nebula (the nominal models of Drouart et al. (1999) and 
Hersant et al. (2001)) intercept the stability curves of the 
different condensates at similar temperature and pressure 
conditions. The stability curve of CO2 pure condensate is 
derived from the existing experimental data (Lide 1999). 
From Fig. ^ it can be seen that CO2 crystallizes as a 
pure condensate prior to being trapped by water to form 
a clathrate hydrate during the cooling of the solar nebula. 
Hence, we assume in this work that solid CO2 is the only 
existing condensed form of CO2 in the solar nebula. 
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mation of Jupiter and its primordial volatile composition, 
"ct^calculated in A05b. These authors showed that, in or- 
i§: ■ to fit the volatile enrichments measured by the Galileo 
HjS-s.TSHjiJ^'be, only some values of C02:CO:CH4 and N2:NH3 
ratios (consistent with ISM observations of Gibb et al. 
' ' "(2(104) and AUamandola et al. (1999)) were allowed. Since 
solids that were incorporated in the Jovian subnebula ini- 
tially formed in Jupiter's feeding zone, they shared the 
same composition than those accreted by proto- Jupiter 
_du:'ing its formation. Hence, in our calculations, we adopt 
'TT^S^HiPB i same C02:CO:CH4 and N2:NH3 ratios in the solar 
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Fig. 1. Stability curves of the species trapped as hydrates 
or clathrate hydrates considered in this work and evolu- 
tionary track of the nebula in P — T space at the heliocen- 
tric distance of 5 AU. Abundances of various elements are 
solar. For CO2, CO and CH4, their abundances are calcu- 
lated assuming C02:CO:CH4 = 30:10:1. For N2 and NH3, 
their abundances are calculated assuming N2:NH3 = 1. 
The condensation curve of CO2 pure condensate (solid 
line) is plotted together with that of the corresponding 
clathrate hydrate (dashed line). The solar nebula cooling 
curve at 5 AU is derived from A05b. 

2.1. Initial ratios of C02:C0:CHi and N2:NH^ in the 
solar nebula gas-phase 

In the present work, the abundances of all elements are 
considered to be solar (Anders & Grevesse 1989) and O, 
C, and N exist only under the form of II2O, CO2, CO, 
CH4, N2, and NII3 in the solar nebula vapor phase. Gas- 
phase abundances relative to II2 in the nebula for species 
considered here are given in Tabled 

Table 1. Gas phase abundances of major species with 
respect to II2 in the solar nebula (from Anders & Grevesse 
for C02:CO:CH4 = 30:10:1 and N2:NH3 = 1. 
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We aim to estimate the composition of ices incorpo- 
rated in Galilean satellites in a way consistent with the for- 



feeding zone 



According to the clathrate hydrate trapping theory 
(Lunine & Stevenson 1985), the complete clathration of 
CO, CH4, N2, NH3, H2S, Xe, Kr, and Ar in Jupiter's 
feeding zone requires an important amount of available 
crystalline water. This tranlates in a H20:Il2 ratio greater 
than that deduced from solar gas-phase abundances of el- 
ements in the solar nebula (see TableQland TableEJ. This 
overabundance may result from the inward drift of icy 
grains (Supulver & Lin 2000), and from local accumula- 
tion of water vapor at radii interior to the water evapora- 
tion/condensation front, as described by Cuzzi & Zahnle 
(2004). The corresponding minimum molar mixing ratio 
of water relative to II2 in the solar nebula gas-phase is 
given by 



XH2O 



E(i?; Th20,Ph2o) r 



(1) 



where Xi is the molar mixing ratio of the volatile i 
with respect to H2 in the solar nebula gas-phase, 7^ 
is the required number of water molecules to form the 
corresponding hydrate or clathrate hydrate (5.75 for a 
type I clathrate hydrate, 5.66 for a type II clathrate 
hydrate, 1 for the NH3-H2O hydrate and for CO2 pure 
condensate), T,{R;Ti, Pi)neb and Y,{R;Th20, PH2o)neb are 
the surface density of the nebula at the distance R from 
the Sun at the epoch of hydratation or clathration of 
the species i and at the epoch of condensation of water, 
respectively. 

Table 12 gives the values of XH2O in Jupiter's feeding 
zone, for the C02:CO:CH4 and N2:NIl3 ratios used in 
A05b and in this work. Note that, in order to calculate 
xh20j we have considered a subsolar abundance for H2S, 
similarly to A05b. Indeed, H2S, at the time of its incorpo- 
ration in icy planetesimals, may have been subsolar in the 
protoplanetary disk, as a result of the coupling between 
the oxygen-dependent sulfur chemistry, the FeS kinetics, 
and the nebular transport processes that affect both oxy- 
gen and sulfur abundances (Pasek et al. 2005). Following 
the calculations described in A05b to fit the observed sul- 
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fur enrichment in Jupiter, we have adopted H2S:H2 = 0.60 
X (S:H2)o for NaiNHg = 10 and HaSiHa = 0.69 x (SiHa)© 
for NaiNHs = 1 in the solar nebula gas-phase. 



The strategy describing the choice of the different sub- 
disk parameters is given in Paper I and the different pa- 
rameters of the Jovian subncbula are recalled in Table |2I 



Table 2. Calculations of the gas-phase abundance of wa- 
ter XH-iO required to trap all volatile species, except CO2 
that condenses as a pure ice, in Jupiter's feeding zone. 
C02:CO:CH4 and NaiNHa gas-phase ratios considered 
here are those determined by A05b in the solar nebula 
gas-phase to fit the enrichments in volatiles in Jupiter. 



NaiNHs = 10 N2:NH3 = 1 

C02:CO:CH4 = 10:10:1 1.55 x lO'^ 1.51 x 10"^ 

C02:CO:CH4 = 20:10:1 - 1.14 X 10"^ 

C02:CO:CH4 = 30:10:1 - 9.48 x 10~* 

C02:CO:CH4 = 40:10:1 - 8.33 x 10"* 



Table 3. Thermodynamical parameters of the Jovian sub- 
nebula. 



Thermodynamical 
parameters 



Mean mol. weight (g/mole) 




2.4 
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X 10-* 


Initial disk's radius 




150 


Initial disk's mass (Mj) 
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X 10~^ 


Initial accretion rate (Mj/yr) 


9 


X 10"^ 



2.3. Composition of ices incorporated in planetesimals 
produced in Jupiter's feeding zone 

Using the aforementioned water abundances, one can cal- 
culate the mass abundances of major volatiles with re- 
spect to H2O in icy planetesimals formed in Jupiter's 
feeding zone. Indeed, the volatile i to water mass ratio 
in these planetesimals is determined by the relation given 
by MG04: 

y- _ ^{R': Tj, Pi)neb 

XH2O TH2O1 PH2o)neb 

where Xi and XH2O are the mass mixing ratios of the 
volatile i and of H2O with respect to H2 in the solar neb- 
ula, respectively. In this calculation, XH2O is derived from 
XH2O in Table 13 

3. Turbulent model of the Jovian subnebula 

3.1. Thermodynamical characteristics of the model 

The a-turbulent model of the Jovian subnebula we con- 
sidered here is the one proposed in Paper I. The subdisk 
evolution is divided into two distinct phases. During the 
first one, the Jovian subnebula is fed by the solar neb- 
ula. During this phase, which lasts about 0.56 Myr, the 
subnebula is in equilibrium since the accretion rate is con- 
stant throughout the subdisk. The origin of time has then 
no influence and is arbitrarily chosen as being the moment 
when Jupiter has already accreted ~ 85% of its total mass. 
When the solar nebula disappears, the accretion rate at 
the outer edge of the subdisk decreases to zero, and the 
subnebula enters its second phase. The subdisk evolves 
due to the accretion of its own material onto the planet, 
and expands outward due to the conservation of angular 
momentum. 



3.2. Evolution of volatile rich planetesimals 
incorporated in the subnebula 

Figure |21 illustrates the fate of ices incorporated in plan- 
etesimals accreted from the nebula by the subnebula. 
From this figure, it can be seen that, as soon as they 
are introduced into the Jovian subnebula, the different 
ices start to vaporize in the whole subdisk. The different 
volatile species considered here start to crystallize again in 
the Jovian subnebula's outer edge between 0.44 and 0.55 
Myr. Note that we have considered here the condensa- 
tion temperatures given in Fig. 1 to calculate the epochs 
of crystallization of the different ices in the subnebula. 
Moreover, we did not take into account the ablation due to 
the friction of planetesimals with gas. Water ice condenses 
at t = 0.57 Myr at the orbit of Callisto (26.6 Rj) and at 
t = 0.61 Myr at the orbit of Ganymede (15.1 Rj). NH3- 
H2O hydrate becomes stable aX t — 0.59 Myr at the orbit 
of Callisto and ait — 0.67 Myr at the orbit of Ganymede. 
CO2 pure condensate is not stable at times earlier than 
0.60 Myr at the orbit of Callisto and 0.68 Myr at the orbit 
of Ganymede. In addition, clathrate hydrates of H2S, Xe, 
CH4, CO, N2, Kr, and Ar become stable between t = 0.59 
Myr and t = 0.68 Myr at the orbit of CaUisto, and t = 0.67 
Myr and t = 0.78 Myr at the orbit of Ganymede. Icy plan- 
etesimals entering into the subnebula at epochs later than 
those indicated above should keep trapped their volatiles 
and maintain the Ices/Rocks (I/R) ratios they acquired 
in the solar nebula. On the other hand, icy planetestimals 
entering into the subnebula at epochs prior to those de- 
termined for preserving ices at the orbits of the two major 
icy satellites must have lost their content in volatiles in 
the satellite zone due to vaporization. 

3.3. Gas-phase chemistry of major C and N bearing 
volatiles in the subnebula 

Since ices were all vaporized in the subdisk during at least 
the first ~ 0.5 Myr of the Jovian subnebula evolution, it 
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Fig. 2. Radii of formation of water ice, NH3-H2O hydrate, 
CO2 pure condensate, and CH4 and Ar clathrate hydrates 
in the Jovian subnebula as a function of time. Radii of 
formation of H2S, Xe, CO, N2 and Kr clathrate hydrates 
are not represented but are within curves 1 and 5. 

seems worthwhile to examine the gas-phase reactions that 
can occur for major C and N volatile species in such an 
environment. 

Following Prinn & Fegley (1989), the net reactions re- 
lating CO, CH4, CO2, N2 and NH3 in a gas dominated by 
H2 are 



CO + H2O = CO2 + H2 
CO + 3H2 = CH4 + H2O 
N2 + 3H2 = 2NH3 



(3) 
(4) 
(5) 



which all proceed to the right with decreasing temperature 
at constant pressure. Reaction Q has been recently stud- 
ied by Talbi & Herbst (2002) who demonstrated that its 
rate coefficient is negligible, even at temperature as high 
as 2000 K (of the order of - 4.2 x IQ-^^ cm^ s'^). Such a 
high temperature range is only reached at distances quite 
close to Jupiter and at early epochs in the Jovian sub- 
nebula (see Fig. 6 in Paper I). As a result, the amount of 
carbon species produced through this reaction is insignif- 
icant during the whole lifetime of the subnebula. 

Reactions Q and Q are respectively illustrated by 
Figs. 13 and ^ The calculations are performed using the 
method described in Mousis et al. (2002a) where the 
reader is referred for details. At the equilibrium, C0:CH4 
and N2:NH3 ratios depend only upon local conditions of 
temperature and pressure (Prinn & Barshay 1977; Lewis 
& Prinn 1980; Smith 1998). C0:CH4 and N2:NH3 ratios 
of 1000, 1, and 0.001 are plotted in Figs.|3|and^ and com- 
pared to our turbulent model at three different epochs (0 
yr, 0.56 Myr and 0.6 Myr). These figures show that, when 
kinetics of chemical reactions are not considered, CH4 and 
NH3 progressively dominate with time in the major part of 
our turbulent model of the Jovian subnebula rather than 
CO and N2. 



However, the actual C0:CH4 and N2:NH3 ratios 
depend on the chemical timescales, which characterize 
the rates of CO to CH4 and N2 to NH3 conversions in our 
model of the Jovian subnebula. We have calculated these 
chemical times from the data given by Prinn & Barshay 
(1977), Lewis & Prinn (1980), and Smith (1998), and 
using the temperature and pressure profiles derived from 
our turbulent model. The results, calculated at several 
epochs of the subnebula's life and at different distances 
to Jupiter, are represented in Fig. El Taking into account 
the kinetics of chemical reactions, one can infer that the 
efficiency of the conversion is limited only to the inner 
part of the Jovian subnebula and at early times of its 
first phase. This implies that C0:CH4 and N2:NH3 ratios 
remain almost constant during the whole lifetime of the 
Jovian subnebula. Moreover, since reaction ij^J plays no 
role in the Jovian subnebula, the C02:C0 ratio also 
remains fixed during its lifetime. 

Finally, these conclusions are compatible with those 
found by MG04 for their colder subnebula model and 
imply that the C02:CO:CH4 and N2:NH3 ratios in the 
subnebula gas-phase were close to the values acquired in 
Jupiter's feeding zone once these species were trapped or 
condensed. From these initial gas-phase conditions in the 
Jovian subnebula, it is now possible to examine the com- 
position of regular satellites ices if these bodies formed 
from planetesimals produced in this environment. 

4. Constraining the composition of ices 
incorporated in regular icy satellites 

Following Paper I, the forming protosatcllitcs migrate in- 
wards the Jovian subnebula under the influence of type I 
migration (Ward 1997). Since protosatellites are suscepti- 
ble to migrate at different epochs of the Jovian subnebula's 
life, two opposite regular satellite formation scenarios can 
then be derived. In the first one, regular icy satellites of 
Jupiter have been accreted from planetesimals that were 
preserved from vaporization when they entered into the 
subnebula after several hundreds of thousands years of ex- 
istence. In the second scenario, regular icy satellites have 
been accreted from icy planetesimals that formed in the 
subnebula. We now explore the consequences of these two 
scenarios on the resulting composition of ices incorporated 
in the Jovian regular icy satellites. 

4.1. First scenario: icy planetesimals produced in the 
solar nebula 

The hypothesis of satellite formation from primordial 
planetesimals (i.e. planetesimals that were produced in 
Jupiter's feeding zone without subsequent vaporization) 
is supported by the recent work of A05b who found, with 
their nominal model for interpreting the enrichments in 
volatiles in Jupiter's atmosphere, that I/R in solids ac- 
creted by the giant planet is similar to that estimated in 
the current Ganymede and Callisto. In this scenario, the 
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Fig. 3. Calculated ratios of C0:CH4 in the Jovian sub- 
nebula at the equilibriimi. The solid line labelled CO-CH4 
corresponds to the case where the abundances of the two 
gases are equal. When moving towards the left side of the 
solid line, CO/CH4 increases, while moving towards the 
right side of the solid line, CO/CH4 decreases. The dot- 
ted contours labelled -3, 0, 3 correspond to logio C0:CH4 
contours. Thermodynamical conditions in our evolution- 
ary turbulent model of the Jovian subdisk are represented 
at three epochs of the subnebula. The Jovianocentric dis- 
tance, in is indicated by arrows when C0:CH4 — 1 
for t = and 0.56 Myr (transition epoch between the two 
phases of the subnebula evolution) . 



mass abundance of major volatilcs with respect to H2O in 
the Jovian regular icy satellites is equal to that in plan- 
etesimals formed in Jupiter's feeding zone, and calculated 
in Sect.ESl 

4.2. Second scenario: icy planetesimals produced in the 
Jovian subnebula 

From Fig. |21 it can be seen that ices entering into the 
Jovian subdisk were all vaporized at epochs prior to ~ 0.5 
Myr. With time, the subnebula cooled down and volatiles 
started to crystallize again following the same condensa- 
tion sequence as that described in the solar nebula (see 
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Fig. 4. Same as Fig.O but for calculated ratios of N2:NH3 
at the equilibrium. The Jovianocentric distance, in is 
indicated by arrows when N2:NH3 = 1 for t = 0, 0.56 Myr 
and 0.6 Myr of our turbulent model. 

Fig.^. One can link the resulting volatile i to water mass 
ratio {Y^s^}j in solids formed into the Jovian subnebula 
to the initial one {Yi)ff,ed in planetesimals produced in 
Jupiter's feeding zone through the following relation: 



(6) 



where fi is the fractionation factor due to the consecutive 
vaporization and condensation of volatile i in the subdisk. 
The fractionation factor fi is given by: 



i Jsub 



sub 



(7) 



where T.{R;Ti, Pi)sub and T,{R;Th20, PH2o)sub are the 
surface densities in the Jovian subnebula, at the distance 
R from Jupiter, and at the epochs of trapping of species 
i and of H2O condensation, respectively. Using conden- 
sation temperatures of the different ices formed in the 
subnebula similar to those calculated in Jupiter's feed- 
ing zone, fi remains almost constant (10 % variations at 
most) in the whole subdisk. Values of fi range between 
0.40 and 0.76 and are given for each species in Table 0] 

In summary, if regular icy satellites were accreted from 
solids produced in the subnebula, the resulting volatile i to 
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Fig. 5. Chemical times profiles calculated for C0:CH4 and 
N2:NH3 conversions in our model of the Jovian subnebula. 
The conversion of CO to CH4 and of N2 to NH3 is fully 
inhibited, except quite close to Jupiter and at the early 
times of the subnebula evolution. 

water mass ratio (Yi)sub hi their ices is that estimated in 
planetesimals formed in Jupiter's feeding zone (see Sect. 
I2.3|l multiplied by the fractionation factor fi. 

Table 4. Mean values of the fractionation factor fi cal- 
culated for icy planetesimals produced in the Jovian sub- 
nebula. 



Species 



fi Species 



NH3:H20 

XeiHaO 

N2:H2 

AriHaO 



0.76 
0.71 
0.59 
0.47 
0.40 



H2S:H20 
CH4:H20 
C0:H20 
Kr:H20 



0.74 
0.57 
0.49 
0.45 



4.3. Composition of regular satellites ices 

Table |S1 summarizes the composition range that can be 
found in the Jovian regular satellites ices formed in the 
framework of the first scenario and assuming that most 
of the trapped volatiles were not lost during the accretion 
and the thermal history of the satellites. From this table, it 
can be seen that C02:H20, C0:H20 and CH4:H20 mass 
ratios vary between 3.7x 10~^ and 1.15, between 1.3x10"^ 
and 1.8 x 10-\ and between 9 x lO'^ and 1.1 x 10"^, 
respectively, in the interiors of regular icy satellites, as 
a function of C02:CO:CH4 and N2:NH3 gas-phase ratios 
assumed in the solar nebula. Similarly, N2:H20, NH3:H20 
and H2S:H20 ratios should be between 3.8 x 10^^ and 
6.9 X 10~^, between 5 x 10~^ and 6.2 x 10~^, and between 
1.9 X 10~^ and 3.6 x 10~^, respectively. Low amounts of Ar, 
Kr, and Xe should also exist in the interiors of regular icy 
satellites. In the second scenario, the resulting volatile i to 
water mass ratio in regular icy satellites must be revised 



down compared to the values quoted above and in Table 
|5l using the fractionation factors given in Tabled 



5. Summary and discussion 

In this work, we have used the evolutionary turbulent 
model of the Jovian subnebula described in Paper I to 
calculate the composition of ices incorporated in the reg- 
ular icy satellites of Jupiter. The model of the Jovian sub- 
nebula we used here evolves in two distinct phases dur- 
ing its lifetime. In the first phase, the Jovian subnebula 
is fed by the solar nebula as long as the latter has not 
been dissipated. In the second phase, the solar nebula has 
disappeared and the subnebula progressively empties by 
accreting its material onto the forming Jupiter. Solids en- 
tering into the Jovian subnebula and that may ultimately 
lead to the Jovian satellite formation are assumed to have 
been produced in the feeding zone of proto-Jupiter prior 
to its appearance. Some of these solids were coupled with 
the material flowing in the Jovian subnebula from the so- 
lar nebula during the first phase of its evolution, due to 
their submeter dimensions, while larger of them, with he- 
liocentric orbits, may have been captured by the subdisk 
when they came through. 

We have considered CO2, CO, CH4, N2, NH3, H2S, 
Ar, Kr, and Xe as the major volatile species existing in 
the gas-phase of Jupiter's feeding zone. All these volatiles, 
except CO2, have been trapped under the form of hydrates 
or clathrate hydrates in Jupiter's feeding zone during the 
cooling of the solar nebula. CO2 crystalhzed as a pure 
condensate prior to be trapped by water and formed the 
only existing condensed form of CO2 in the feeding zone 
of Jupiter. 

We employed C02:CO:CH4 and N2:NH3 ratios consis- 
tent with those used by A05b, namely C02:CO:CIl4 be- 
tween 10:10:1 and 40:10:1, and N2:NH3 between 1 and 10 
in the gas-phase of Jupiter's feeding zone. Such a range 
of values is compatible with those measured in ISM or 
estimated in the solar nebula. This allowed us to deter- 
mine the corresponding minimum Il20:Il2 gas-phase ra- 
tios required to trap all volatiles (except CO2) in the giant 
planet's feeding zone. 

Moreover, since, according to our model, ices contained 
in solids entering into the subnebula before ^ 0.5 Myr 
were all vaporized, we have followed the net gas-phase 
chemical reactions relating CO, CII4, CO2, N2, and NII3 
in this environment. We then concluded that these reac- 
tions are mostly inefficient in the Jovian subnebula, in 
agreement with the previous work of MG04. This involves 
that C02:CO:CIl4 and N2:NIl3 ratios were not essen- 
tially different from those acquired in the feeding zone 
of Jupiter, once these species were trapped or condensed 
in the subdisk. In addition, in order to estimate the mass 
abundances of the major volatile species with respect to 
II2O in the interiors of the Jovian regular icy satellites, 
we considered the formation of these bodies by following 
two opposite scenarios. 
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Table 5. Calculations of the ratios of trapped masses of volatiles to the mass of H2O ice in regular icy satellites 
accreted from planetesimals formed in Jupiter's feeding zone. Gas-phase abundances of H2O are given in Table 2 and 

gas-phase abundances of elements, except S (sec text), arc assumed to be solar (sec Table 1). Ranges of C02:CO:CH4 
and N2:NH3 gas- phase ratios considered here are those determined by A05b in the solar nebula gas-phase to fit the 
enrichments in volatiles in Jupiter (see text). 



N2:NH3 = 10 NaiNHg = 1 



Species 


C02:C0 = 


: 1 


C02:C0 = 


: 1 


C02:C0 = 


: 2 


C02:C0 = 


3 


C02:C0 = 


: 4 


C02:H20 


3.72 X 10" 


-1 


3.84 X 10" 


-1 


6.92 X 10" 


-1 


9.43 X 10" 


-1 


1.15 




C0:H20 


1.78 X 10" 


-1 


1.83 X 10" 


-1 


1.63 X 10" 


-1 


1.47 X 10" 


1 


1.34 X 10" 


-1 


CH4:H20 


1.14 X 10" 


-2 


1.17 X 10" 


-2 


1.04 X 10" 


-2 


9.44 X 10" 


■3 


8.60 X 10" 


-3 


NaiHaO 


5.34 X 10" 


-2 


3.82 X 10" 


2 


5.06 X 10" 


-2 


6.07 X 10" 


-2 


6.90 X 10" 


-2 


NH3:H2 


4.60 X 10" 


■3 


3.43 X 10" 


-2 


4.55 X 10" 


-2 


5.45 X 10" 


-2 


6.20 X 10" 


-2 


H2S:H20 


1.93 X 10" 


2 


1.99 X 10" 


-2 


2.63 X 10" 


-2 


3.16 X 10" 


2 


3.59 X 10" 


-2 


Ar:H20 


4.17 X 10" 


-3 


4.29 X 10" 


-3 


5.69 X 10" 


3 


6.83 X 10" 


-3 


7.76 X 10" 


-3 


KriHaO 


4.89 X 10" 


-6 


5.04 X 10" 


-6 


6.68 X 10" 


■6 


8.01 X 10" 


■6 


9.11 X 10" 


-6 


XerHaO 


9.29 X 10" 


-7 


9.57 X 10" 


■7 


1.27 X 10" 


-6 


1.52 X 10" 


-6 


1.73 X 10" 


-6 



In the first scenario, regular icy satellites were ac- 
creted from planetesimals that have been preserved from 
vaporization during their migration in the Jovian subncb- 
ula. This assumption is in agreement with the work of 
A05b who found, with their nominal model for interpret- 
ing the enrichments in volatiles in Jupiter's atmosphere 
(N2:NH3 = 1 and C02:CO:CH4 = 30:10:1 in the solar 
nebula gas-phase), that I/R in planetesimals accreted by 
the giant planet is similar to those estimated by Sohl ct 
al. (2002) in Ganymede and Callisto. This allowed us to 
estimate the ratios of the trapped masses of volatiles to 
the mass of H2O ice in the regular icy satellites, assuming 
that these species were not lost during their accretion and 
their thermal history. 

In the second scenario, regular icy satellites were ac- 
creted from planetesimals produced in the subnebula. 
Indeed, in the framework of our model, ices contained 
in solids were entirely vaporized if they entered at early 
epochs into the Jovian subdisk. With time, the subnebula 
cooled down and ices crystallized again in the subneb- 
ula prior to having been subsequently incorporated into 
the growing planetesimals. In this second scenario, assum- 
ing, as in the first one, that the regular icy satellites did 
not lose volatiles during their accretion phase and their 
thermal history, we have also estimated the composition 
range of ices trapped in their interiors. In that scenario, 
the amount of icc^s incorporated in regular icy satellites 
should be lower than in the previous one where planetes- 
imals were produced in the solar nebula. 

In both scenarios, the calculated composition of the 
Jovian regular satellites ices is consistent with some ev- 
idences of carbon and nitrogen volatile species in these 
bodies, even if the presence of some predicted compo- 
nents has yet to be verified. For example, reflectance spec- 
tra returned by the Near-Infrared Mapping Spectrometer 
(NIMS) aboard the Galileo spacecraft revealed the pres- 
ence of CO2 over most of Callisto's surface (Hibbitts et 
al. 2000). Moreover, Hibbitts et al. (2002) suggested that 



CO2 would be contained in clathrate hydrates located in 
the subsurface of Callisto and would be stable over ge- 
ologic time unless exposed to the surface. CO2 has also 
been detected on the surface of Ganymede (Hibbitts et 
al. 2003). In addition, one explanation for the internal 
magnetic fields discovered in both Ganymede and Callisto 
(Kivelson et al. 1997, 1999) invokes the presence of sub- 
surface oceans within these satellites (Sohl et al. 2002). 
The presence of such deep oceans is probably linked to 
the presence of NH3 in the interiors of these satellites, 
since this component decreases the solidus temperature 
by several tens of degrees (Mousis et al. 2002b; Spohn & 
Schubert 2003). 

Subsequent observations are required to determine 
which of both presented formation scenarios is the most 
realistic. On the basis of isotopic exchanges calculations 
between HDQ and H2 in the solar nebula, Mousis (2004) 
estimated that the D:H ratio in the Jovian regular icy 
satellites is between ~ 4 and 5 times the solar value, as- 
suming they were formed from planetesimals produced in 
Jupiter's feeding zone. On the other hand, icy satellites 
formed from planetesimals produced in the Jovian sub- 
nebula should present a lower D:H ratio in H2O ice since 
an additional isotopic exchange occurred between HDO 
and H2 in the subdisk gas-phase. Such estimates, com- 
pared with further in situ measurements of the D:H ratio 
in H2O ice on the surfaces of the Jovian regular satellites, 
should allow to check the validity of the two proposed 
formation scenarios. 
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